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^ ■ ABSTRACT 

o , 

CO , Non-parametric, empirically based, models for associating galaxy luminosities 

with halo/subhalo masses are being developed by several groups and we present here 
an updated version of the lVale fc Ostrikeil l|2004f) version of this model. This is based 
on a more accurate, self-consistent treatment of subhalo mass loss and revised results 
for the subhalo mass function to address this question anew. We find that the mass- 
QQ \ luminosity relation, at high mass, particularly for first brightest galaxies and less so 

for group total, is almost independent of the actual luminosity function considered, 
when luminosity is scaled by the characteristic luminosity L*. Additionally, the shape 
i of the total luminosity depends on the slope of the subhalo mass function. For these 

■ high mass, cluster sized haloes we find that total luminosity scales as Ltot ~ M°-^*, 

while the luminosity of the first brightest galaxy has a much weaker dependence on 
Q^. halo mass, Li M°-^^ , in good agreement with observations and previous results. 

1 ' At low mass, the resulting slope of the mass luminosity relation depends strongly of 
O I the faint end slope of the luminosity function, and we obtain a steep relation, with 

approximately L ~ M'*-^ for AI ^ lO^°/i~^M0 in the K-band. The average number of 
^ ' galaxies per halo/cluster is also in very good agreement with observations, scaling as 

In general, we obtain a good agreement with several independent sets of obser- 
vational data. Taking the model as essentially correct, we consider two additional 
possible sources for remaining discrepancies: problems with the underlying cosmology 
' and with the observational mass determination. We find that, when comparing with 

observations and for a flat cosmology, the model tends to prefer lower values for fi^ 
and as- Within the WMAP-I-SDSS concordance plane of iTegmark et alJ l)2004D . we 
find best agreement around ilm — 0.25 and erg = 0.8; this is also in v ery good agree- 
ment with the results of the CMB-|-2dF study of lSanchez et alJ l)2005j) . We also check 
on possible corrections for observed mass based on a comparison of the equivalent 
number of haloes/clusters. Additionally, we include further checks on the model re- 
sults based on the mass to light ratio, the occupation number, the group luminosity 
function and the multiplicity function. 

Key words: galaxies: haloes - cosmology: theory - dark matter - large-scale structure 
of the universe 



1 INTRODUCTION 

One of the outstanding challenges of cosmology is to re- 
late the more theoretical aspect of the standard cosmologi- 
cal model, in the form of the large scale distribution of dark 
matter as seen in high resolution N-body simulations, to 
the observational evidence, as reflected by large galaxy sur- 
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veys. Or, in a directly related question, what is the relation 
between observable properties of galaxies and computable 
properties of dark matter haloes? 



The traditional route to looking at this problem has 
been to follow the theory of galaxy formation. This can be 
tested through the results of hydrodynamical simulations 
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I Someryille et alJ l200lt IWechsler et alJ I2001I : iBenson et all 
l20n,'j:lBerlind et al.ll20n.'j^ . This is a powerful way of looking 
at the problem, since it provides a direct answer, and a 
galaxy formation theory is an end goal in itself. There are, 
however, some problems: first and foremost, the theory 
behind galaxy formation has several components which 
are ill understood, and where the best approach remains 
phenomenological. At the same time, the required very high 
resolution, large scale, full hydrodynamical simulations 
including magnetic fields and radiative transfer are beyond 
current computing resources. 

In the past few years, a new approach has appeared 
which is in many ways an alternative: this consists of 
an indirect method, associating galaxy and dark matter 
halo properties in an empirical manner. This is usually 
done through a statistical approach, either by focus- 
ing directly on the number of galaxies in each h alo, as 
is do ne in halo occup ati on distribution models iSeliak 
2OOOI : iBensonI I2001I: iBullock. Wechsler. fc Somervilk 
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or through the luminosity distribution of the galax- 
ies in a halo, like the conditional luminosi ty func- 
tion approach llvan den Bosch. Yan g fc Mol 120031: 
[^anEi^^^^^^ai^den"Boscl7~ '2003^: I van den Bosch et alJ 
I2OO5I: ICoorav fc M ilosavlicvic 2005bi ICooravl l2005dlb^ . or 
finally by building a direct relation between mass and 
luminosity JPeacock fc SmithI I2OO0I: iKravtsov et alJ 12004 
I Vale fc Qstrike jl2004l : iTasitsiomi et alJl2004) . 

In the present paper, we take the latter approach. We 
use new, high resolution simulation results for the subhalo 
mass distribution, together with a self consistent approx- 
imation to subhalo mass loss, to build the total distribu- 
tion of dark matter hosts based on physical theory. To- 
gether with the empirically determined galaxy luminosity 
function, this yields a mass luminosity relation for an indi- 
vidual galaxy and the halo or subhalo which hosts it. This is 
in part an updated versio n of the preliminary model shown 
in I Vale fc Ostrikeij (|2004 hereafter paper I). 

This paper is organised as follows: in section 2, we in- 
troduce the non-parametric model and the main concepts 
behind it. In section 3, we look in some detail at the subhalo 
mass function, introducing a prescription for subhalo mass 
loss and in particular looking at how to define this mass 
function in terms of the original, pre-merger into parent, 
mass of the subhaloes. In section 4, we give the main results 
for our base model, using the K-band luminosity function 
from the 2MASS survey. In section 5, we look at how the 
model results change depending on the luminosity function 
(and hence waveband) used, as well as the effect of chang- 
ing the underlying cosmological model and what we find is 
the most relevant parameter of the subhalo mass function, 
the low mass slope. We also look at potential problems with 
the observational mass determination. We present additional 
checks of the model in section 6, by calculating the predicted 
mass to light ratio, occupation number, group luminosity 
function and multiplicity function. Finally, we conclude in 
section 7. 



2 THE NON-PARAMETRIC MODEL 

The basic idea behind the non-parametric model is to as- 
sume a monotonic and one to one relation between the lu- 
minosity of a galaxy and the mass of the dark matter halo 
or subhalo which hosts it. This comes from the standard 
picture that galaxies are formed in haloes through the ac- 
cretion of gas, the amount of which is a monotonic function 
of the depth of the potential well of the halo and thus of 
its mass. By including the subhaloes, we are assuming that 
all galaxies are either hosted individually in a parent halo, 
or in the case of multiple systems like clusters in one of 
the subhaloes. In fact, we assume that counting haloes and 
subhaloes accounts for all hosts, and exclude complications 
such as conditions during ongoing mergers. Groups and clus- 
ters are then formed when haloes merge; the end result of 
such a build up is to have a central galaxy, which formed in 
the most massive of the initial haloes (which subsequently 
became the parent), and satellite galaxies in the smallest 
haloes which were accreted and which became subhaloes. 

Using this basic concept, it becomes possible to ob- 
tain an average relation between the luminosity of a galaxy 
and its host halo/subhalo mass by matching the numbers of 
each, the former from large scale galaxy surveys, the latter 
from dark matter simulations. The average luminosity L of 
a galaxy in a halo or subhalo of mass A4 will then simply be 
given by: 



4>{L)dL= I n(M)dM , 



(1) 



where 4'{L) is the galaxy luminosity function and n{M) = 
nn^M) + nsH{M) is the sum of the halo and subhalo mass 
functions. 

During and after the halo merging pro- 
cess significant star formation in the subh al o de- 
clines (e.g., |Kaufl hnajarij_WIiit^^ Il993l : 
ISomerville fc Primackl Il999l: ICole et al.l bOOolT so we 

will take as a constraint that the amount of gas accreted by 
what becomes a subhalo, and which will thus be available 
to form a satellite galaxy, will be proportional not to the 
mass it has at present, but to the maximum mass it had 
before being accreted by its parent halo and undergoing 
mass loss to tidal stripping and dy namical friction. This is 
supported by some recent work of 'Libesk ind et all i2005l) . 
who use high resolution N-body simulations together with 
semi-analytical modelling of the formation of galaxies to 
study the distribution of satellites in Milky Way type 
galaxies. They find that while the spatial distribution of 
satellite galaxies is significantly different from that of the 
most massive present day subhaloes, it is well matched 
by that of the subset of subhaloes with the most massive 
progenitors. It is then necessary to use not the distribution 
of subhaloes as a function of their present mass, but instead 
of the original mass they had prior to accretion into the 
parent and subsequent mass loss, regardless of theoretical 
issues related to galaxy formation. We circumvent this 
problem by coupling the present, evolved subhalo mass 
function with a prescription for the amount of mass loss, 
together with some simple arguments on the total mass 
contained in these subhalo progenitors, in order to regain 
the initial mass of a subhalo at the moment of accretion 
and tag each of these by its initial mass. 



The Non-Parametric Model 3 




0.0001 



Figure 1. Mass loss factor f {m„ow / 'mo) y which is a function 
both of subhalo and parent halo mass. Plotted are curves for 
three different parent halo masses, 10^^/i~^Mq, 10^^h~^MQ, 
10^^ h~^MQ, from top to bottom respectively. These curves have 
th e normalization we obtained direct l y from adopting the results 
of lvan den Bosch. Tormen &: Giocolil 12005 ^: it may be necessary 
to renormalize them when considering different present subhalo 
mass functions (see text for discussion). 



3 SUBHALOES 
3.1 Mass loss 

As discussed, the first step to having a workable non- 
parametric model is to account for the amount of mass 
loss in subhaloes. To obtain it, we ad apt the results of 
Ivan den Bosch. Tormen fc Giocolil ll2005ll , by comparing the 
results they give for the evolved and unevolved (original) 
subhalo mass functions. If we assume that the present and 
original masses of the subhaloes are related monotonically 
(that is, that the most massive of the original subhaloes are 
still the most massive at present, independently of the ac- 
tual amount of mass loss), we can obtain this relation by 
comparing the total numbers of present subhaloes to their 
original progenitors, in a manner similar to the process pre- 
sented above for the non-parametric model. 

The resulting mass loss factors will depend on both 
the subhalo mass and the mass of the parent halo, as can 
be seen in figure Q In fact, we have made a small further 
change to the results wc obtain from the mass functions of 
[van den Bosch. Tormen & Giocoli (2005), which is to flat- 
ten the mass loss factor to a constant at its minimum and 
so maintain monotonicity. The upturn at high subhalo mass 
is explained by two competing factors which contribute to 
mass loss: while the more massive subhaloes lose a larger 
fraction of their mass in each orbit within the parent halo, 
they will also, on average, have formed and therefore been 
accreted later, which means they will have undergone fewer 
of these orbits than less massive subhalo es (see for example, 
Ivan den Bosch. Tormen fc GiocolilboOSh . 



We can now apply this mass loss factor to subhalo mass 
functions measured from simulation results to obtain the 
original subhalo mass function. There is an important fur- 
ther consideration, though: care must be taken with the nor- 
malization of this function. In fact, we know that the sum 
of the mass in these original subhaloes must equal the total 
present mass in the parent halo, since the latter is built up 
by accreting and stripping subhaloes (ignoring for simplic- 
ity the small amount of mass associated with subhaloes that 
disappear totally). Therefore, when using different subhalo 
mass functions, the result presented in figure should be 
taken only to be the shape of the mass loss factor. It is then 
necessary to either renormalize the mass loss factor (or alter- 
natively the derived original subhalo mass function) to en- 
sure that the total mass in these original subhaloes matches 
the present mass of the parent halo. 



3.2 The subhalo mass function 

The distribution of the subhaloes is based on the subhalo 
mass function (SHMF), N{m\M)dm, which gives the num- 
ber of subhaloes in the mass range m to m -I- dm for a parent 
halo of mass A/. We start by defining the present day mass 
fraction in subhaloes as: 



7(M) . - 



mN{m\M)dm, 



(2) 



where mcut{M) is the mass of the most massive subhalo 
possible, which we take to be a function of the parent halo 
mass. This value will be set by assuming that the maximum 
original mass is half the total mass of the parent, M/2. This 
is in fact a question of definition, since a subhalo which had 
a mass greater than this would actually be larger than any 
other it could be merging with, and would itself in fact be 
the parent halo. This can be converted back to rricut by using 
the mass loss factor. 

This mass fraction in subhaloes is generally found to 
be a growing function of parent halo m ass, both in sim- 
ulations or semi-analytical models (e.g .j_|Ga^^^^Jj2004; 
Ivan den Bosch. Tormen fc Giocoli"2005: 'Shaw et al.'*200^. 
As pointed out in van den Bosch, Tormen fc Giocoli (2003), 
this is naturally due to the later formation of more massive 
haloes, therefore allowing less time for mass loss to occur. 
There also seems to be some agreement that the value of 
this mass fraction should be just under 10%, even at high 
halo mass. This will have as an important consequence that 
the subhalo mass function cannot possibly be universal but 
must depend on M, since the subhalo mass fraction in effect 
sets its normalization. 

Following wh at was found by IWeller et alJ i2004l) : 
IShawet ai]ll2005l) . we can then take the SHMF to have the 
form of a Schechter function: 



N{m\M)dm = Am 



Am = 



7(M) 



exp 



m \ dm 



P{V[2 - q] - r[2 - a, m,^t[M)/pM]) 



(3) 



(4) 



where a gives the low mass slope and [3 represents 
an additional cutoff mass. Typica l values for the low 
mass slop e a are around ~ 1.9 Lucia et alJ 
iGao et al.l 12004 | van den Bosch. Tormen fc Giocolil 
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2005 
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IZentner et al.l2005l) . althoueh lshaw et alj ll2005l) have found 
a less steep a = 1.75 when fitting a power law and an 
even flatter a = 1.5 when fitting to a Schechter func- 
tion of the form of equation ijH. We should a lso point 
out that Ivan den Bosch. Tormen fc Giocoli ll2005h obtain a 
slope which varies slightly with parent halo mass. The cutoff 
value has alternatively been determined as = 0.13 from 
semi-analytical models jvan den Bos ch. Tormen fc Giocoli 
[ 2005 ), or /3 = 0.3 from simulations iShaw et ah 2005). To 
this expression should be added the cutoff we are implicitly 
considering at rricut, as expressed in equation The first 
term on the right hand side guarantees that the mass frac- 
tion in subhaloes is 7(M). From this we see that rricut has 
essentially an effect on the normalization of the SHMF. 

However, as discussed above, what we really require 
to build the non-parametric model is the original subhalo 
mass function instead. To do this, we apply the mass loss 
factor presented in the previous section to the SIfMF of 
equation 0, that is, we use m„ow = f{m„om, M)mo, with 
firrinow, M) the mass loss function, renormalizing as appro- 
priate to guarantee that the original total mass in subhaloes 
equals the present mass of the parent halo. In terms of the 
original mass mo, the SHMF of equation |3 then becomes 



No{mo\M)dmo = N[f{m, M)mo\M]f{m, M) x 

dlog/KA^)\ 

Ji ""^0 ' 

dlogm / 



(5) 



where /(m, M) is the mass loss factor, with m the present 
subhalo (i.e., the function plotted in figure^; we have also 
explicitly included the transform of the differential term dm. 
Once more, we are implicitly assuming a cutoff at an original 
subhalo mass of m = M/2. FigureHshows different SHMFs 
calculated using this scheme. 

There are some important conclusions to be drawn from 
figure 13 First, unlike the case with the present SHMFs, 
the original ones are almost universal when plotted as a 
function of m/M. This is a result of the fact that we are 
putting in what is, in these units, a universal normaliza- 
tion at the total halo mass M and also a universal cutoff at 
m/M = 0.5, which usually dominates over the one in equa- 
tion after transforming to the original mass. As a direct 
consequence of this, we can see that the original SHMF, 
as calculated through this method, is essentially indepen- 
dent of the present subhalo mass fraction 7 and largely in- 
dependent of the cutoff given by the parameter /3. It is also 
largely insensitive to the shape of the mass loss factor as 
determined in the previous section. However, it should be 
cautioned that this depends on the actual mass loss factor 
/ used: this is only true if it is fairly regular, or more pre- 
cisely dlog//dlogm ^ 1, as is mostly the case of the function 
we are using presently, which mostly causes a rescalling of 
the subhalo mass function, without much altering its shape. 
In such a case, in fact, the only parameter of the present 
SHMF which has a large effect is the low mass slope a. 
This indicates that it should be mostly correct, in the con- 
text of building a non-parametric model like we are doing 
here, to use as the original SHMF a power law of slope a, 
with a cutoff at m/M = 0.5 and the whole normalized so 
that the total original mass in subhaloes equals the present 
parent halo mass. From this, it is important to retain that. 
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Figure 2. The subhalo mass function (SHMF): shown are 
the original SHMFs for parent haloes of mass IQi^h-lM© and 
lO^^/i~^M0. For comparison, we also show the original SHMF for 
a IQl^h-lMo halo when we take a flat mass loss factor, instead 
of one based on the curves of figure and a present, evolved 
SHMF with a subhalo mass fraction of 7 = 0.1. 



accepting these two conventions we are using on the normal- 
ization and cutoff, the only element of the SHMF to which 
our model is sensitive is the low mass slope. 

For the subsequent calculations, we will use a model 
with the mass loss factor calculated in the previous section, 
a present SHMF given by equation @ with a = 1.9 and 
f3 — 0.3. As mentioned above, the actual value of the mass 
fraction 7 is not important since we have to renormalize the 
original mass function we obtain; it is only factored if we 
wish to have the appropriately normalized mass loss factor 
fim„ow, M) to convert between original and present subhalo 
mass using mnow = fijnnow, M)mo. Hereinafter, unless oth- 
erwise stated, all subhalo masses refer to the original mass. 



4 BASIC MODEL 

In this section we present our base model. We take a flat 
cosmology with parameters Oa/ = 0.25, erg = 0.8 and 
h — 0.7. While these do not exactly corr e spond to the 
current standard model ([Bahcall et al.l [19991 ISoergel et alJ 
l2003t iTeemark et al.l l2004l) . thev are within the aUowed 
range of the flm — gg plane determined fr om the joint 
WMAP-SDSS study of iTegmark et alj (|200i), and we find 
they produce results better matching observations (see sec- 
tion l5.2l below) than do models with slightly higher {Qm, o"g). 
On the other hand, they match very well with the results 
found from an analy sis of CMB and 2dF power spectrum by 
ISanchez et al.l 1I2OO3), being near the center of their — as 
concordance region (see figure]^. For the basic model, we 
use the luminosity function in the K-band a s determined 
from the Two Micron All-Sky Survey f2MASS: Uarrett et alJ 
I^OOs). Using the K-band allows us to avoid possible com- 
plications and discrepancies arising from brief intervals of 



The Non-Parametric Model 5 



1.2 

1.1 

1.0 
0.9 



« 0.8 

b 



0.7 



0.6 



0.5 



0.4 




WMAP+SDSS Tegmaricetal. 2004 

a+^dF Sanchez etal, J(JU5 
Hoekstriiet.il 2002 
Base model 



0.1 0.2 0.3 0.4 0.5 0.6 

Figure 3. The Qm — piano. Shown are the 1- and 
2(T countours from the CMB+2dF power spectrum results of 
LSanchcz ct al. 12005), together wit h the 2cr concordance region 
from the WMAP+SDSS results of iTegmark et alj <2004 . Also 
sh own as an example of weak lensing results are the constraints 
of iHoekstra. Yee. &: GladdersI l2002h . The point marks the pa- 
rameter values chosen for our basic model. 

active star formation, and therefore helps present a clearer 
picture. We use here a single waveband for clarity in the con- 
struction of the model; a comparison of results for different 
wavebands using different luminosity functions is presented 
in section I^TTI 

4.1 Host mass function 

The first step to building the non-parametric model is to 
combine the above subhalo mass function with the parent 
halo distribution in the form of the halo mass function to 
obtain the global distribution of the dark matter hosts of 
galaxies. This is also the step where the cosmological pa- 
rameter dependency is factored in, by determining the form 
of the mass function. We will use the Sheth-Tormen mass 
function (Shcth & Tormen 1999), given by: 



nh{M)dM = A 



with u 



hi) 



2 du 
nWdM 



exp 



dM, (6) 



■'D(.)aiM)^ « = 0-™7, A ^ 0.322 and q = 0.3; 
as usual, o{M) is the variance on the mass scale M, D{z) is 
the growth factor, and 5c is the linear threshold for spherical 
collapse, which in the case of a flat universe is 5c = 1.686. 
This mass function can be roughly approximated by a power 
law at low mass, scaling as ~ M~^'^*, and an exponential 
cutoff at high mass of the form exp{—M/M,), where the cut- 
off mass is defined by M/Mt, = /2 and we have roughly 
M* 5 X lO"/i~^M0. The total distribution of dark mat- 
ter hosts is then given by combining this with the global 
distribution of subhaloes, which can be calculated from the 
SHMF ^ and the halo mass function Q by 
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Figure 4. Global mass functions for haloes and subhaloes. In the 
case of the subhaloes, the x-axis refers to the original, unstripped 
mass (see discussion in text). 
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nh{M)N{m\M)dM . 



(7) 



These mass functions are shown in figure |1] Note that in 
the combined distribution of haloes/subhaloes the former 
dominate with subhaloes approaching parity only on mass 
scales below 10"/i~^Mn. 



4.2 Luminosity Function 

The galaxy distribution is accounted for by the galaxy lu- 
minosity function, which wc assume takes the shape of a 
Schechter function (Schcchtcr 1976): 

^(L)dL = 0.(i-)%xp(-il)^. (8) 

For this base model, we use the K-band results from the 
2MASS survey, with parameters given by a = —1.09, 
= 1.16 X 10~^fc^ Mpc~^ and M, - Slog/i = -23.39 
Kochanek et alj|200 j) (see also table 1). In section 15. II be- 
low, we analyse the result of using luminosity functions in 
other wavebands. This fit applies to the magnitude range 
—22 > Mk — blog^gh > —25. An important point to note is 
that we will assume that the fit actually holds outside of this 
limited range, and extrapolate its result outside of it as nec- 
essary. This raises the question of how good a description of 
the real luminosity function a Schechter function is and, in 
particular, whether a power law is a good fit to the faint end. 
Observations at the faint end are plagued by a host of natu- 
ral difficulties, most of all the need to sample efficiently low 
luminosity and low surface brightness galaxies. It is not yet 
entirely certain whether the luminosity function at the faint 
end can be well fit by the power law in the Sch echter fu nction 
over a wide range (see for example [ Trentham & Tull^|200^ 
[TYentham. Sampson fc Baneriill2005l : see also iBlanton et alJ 
I2OO5I who look at faint galaxies in SDSS and find that a 
Schechter function is a poor fit, and need to introduce a 
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Figure 5. Mass K-band luminosity relation, both for single galax- 
ies and total group luminosity (upper curve). Galaxy luminos- 
ity applies to both haloes and subhaloes, where the mass of the 
latter refers to the original, prior to a ccretion into the p arent 
case. The data points s hown come from lLin fc Motul ^2004^ and 
iLin. Mohr fc Stanford! i2004h (on cluster scales, for both total 
and b rightest galaxy luminosity), analysis of SDSS weak lensing 
iYang et al. 2003; Coorav & Milosavlicvic 2005a), and an esti- 
mate of Milky Way dwarf spheroidal mass and luminosity (see 
text for details). 



double rather than single power law parametrization at the 
faint end). For now, we will simply assume that we can ex- 
tend the measured faint-end slope of the 2MASS survey to 
fainter magnitudes than their stated limit. 



4.3 Mass Luminosity relation 

We can now combine the total dark matter host distri- 
bution with the galaxy distribution to get a relation be- 
tween halo/subhalo mass and galaxy luminosity, by match- 
ing counts as described above in section 2. The resulting 
relation is shown in figure |S| 

Also shown are observational data for brightest clus- 
ter galaxy luminosity and cluster lumi nosity versus cluster 
mas s, taken from t he K-band results of iLin fc Mohrl ||2004|) 
and lLin. Mohr fc~S tanford (2004). These studies use X-ray 
temperature results to determine the cluster mass together 
with luminosities taken from the 2MASS survey; we include 
a correction to the halo mass to account f or a virial overden- 
sity o f 100 times the critical density (e.g.. lBrvan fc NormanI 
Il998l) . instead of 200, based on taki ng an NFW halo with 
a conc entration of 5 as was done in lLin. Mohr fc StanfordI 
(|20o3) (see further discussion in section 15.31 . The points 
at the low-end come from Milky Way dwarf spheroidal 
satellites, where the mass was estimated from results of 
[Havashi e t al. ( 200j) (see also discussion in p aper I) and th e 
luminosity comes from the B-band results of iMateol il99dl . 
then using an average B — K = 3 .6 to convert to the K- 
band iMobasher &: Trenthamll998l) . The intermediate mass 
results come from an analysis of SDSS weak lensing results 



in lYang et all (l2003l) and ICoorav fc Milosavlie^ (l2005bl) ; 
in both of these, we have transformed the z-band luminosity 
to K-band by using log(Lx/ fz) = -0.014 -I- 0.492(g - r), 
with an average {g-r) = 0.6 iBell et alj|2003l) . For all these 
cases, w e again conver t ed to Mioo, using for the concentra- 
tion the lBullock et alJ ll200j) model. 

We find that the mass luminosity relation we obtain 
(figure!^ can be well approximated by a double power law 
of the type: 



{m/m'Y 



(9) 



with the parameters Lo = 1.23x10^°, m' = 3.7x lO^/i'^M©, 
a = 29.78, b = 29.5, and k = 0.0255. This fit is good to 
within ~ 6% in the range 5 X 10^" < M/{h-^MQ) < 5x10^^ 
The most important fact to retain from this fit is that lu- 
minosity scales as L oc M°-2» at high mass. The high values 
obtained for the exponents a and b are an artifact of the 
fact that the relation is steepening as the mass decreases 
(note the value of the break mass given by m'). Group lu- 
minosity scales SlS -Li group CXl M°-^^ for high halo mass. Over- 
all, these parameters look different from the fit to observed 
da ta done using the same functio nal form as equation Q 
bv lCoorav fc Milosavlie^ l^2005b^ . The main reason for the 
difference is the corrections we put in for the mass estimates. 
Also, as noted, care should be taken with the low mass slope, 
since the values for the best fit curve depend on the lowest 
mass used for the fit. It is, however, easy to predict this 
slope: if the mass function goes as ~ M~°'dM (where this 
represents the total for haloes plus subhaloes) and the lumi- 
nosity function as ~ L^^dL, then the luminosity will scale 
as L cx M", with c = (a - l)/(6 - 1). With a ~ 1.95 and 
b = 1.09, we get c ~ 10.5. 

Overall, these results seem a good match to the obser- 
vations. The quality of this match has improved over the re- 
sults we presented in paper I, especially for the group/cluster 
luminosity. As we had hinted in the discussion in paper I, 
this is most likely due to the fact that we now have a bet- 
ter, more self consistent, way of treating the subhalo mass 
fraction and mass loss. It is also interesting to note that, 
since we obtain a total luminosity which scales almost lin- 
early with mass, we naturally obtain a mass to light ratio 
on cluster scales that is almost constant; this is intrinsic to 
the model, and is in good agreement with observations (see 
section 16.11 below) . It is also worth to point out that the 
agreement for the group luminosity is not trivial: only the 
halo/subhalo mass-galaxy luminosity relation is constructed 
directly from the model. The group luminosity is obtained 
by assuming that applying this derived galaxy luminosity 
to the system of the haloes and their subhaloes will result 
in a good match for the corresponding galaxy systems. The 
agreement we obtain seems to point out that this is indeed a 
valid assumption (but see Coorav & Con 2005 for some po- 
tential problems involving the detailed luminosity functions 
for individual groups when following this prescription). 

There is some discrepancy between our results and 
the results for the higher mass bins based on weak lens- 
ing of SDSS galaxies (the po ints from [Yang et al. 2 00^ : 
ICoorav fc Milosavlievidl2005al) . While there is no straight- 
forward way to explain the differences we find, we caution 
that there are several factors that, when taken together, may 
be enough to account for them. Most relevant of these are 
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the fact that we are using an average colour to change the 
luminosities from the z- to the K-band, and the possibility 
that the mass being accounted for does not correspond ex- 
actly to the definition we are using. Intriguingly, these weak 
lensing results seem to be a good match to our curve for 
the total luminosity, which raises the additional possibility 
that the observed luminosity includes not only the central 
lensing galaxy, but also some close, faint, satellites (since we 
would expect that most substructure in these medium mass 
parents to be relatively small). 



5 DEPENDENCE ON LF AND COSMOLOGY 
5.1 Luminosity function 

In this section, we look at how our base result for the mass 
luminosity relation changes depending on the luminosity 
function (and consequently waveband) we use as a starting 
point. We use a variety of published luminosity functions, 
listed in table 1, ranging from the blue to the infrared, and 
repeat the analysis of the previous section for each of these. 
The results for the galaxy and total luminosity are shown 
in figure |H| where the luminosity is plotted in units of the 
characteristic luminosity L* of each of the luminosity func- 
tions. 

The most striking feature of these results is that galaxy 
luminosity in high mass haloes is almost independent of the 
luminosity function used, with only very slight differences 
due to different normalization and faint end slope. Techni- 
cally, this is a consequence of the fact that all luminosity 
functions considered have a bright end cutoff of the form 
exp{—L/Lt,). In general, galaxies have an observed colour- 
magnitude relation which should be reflected in a different 
mass-luminosity relation at different wavebands (as is the 
case here for lower mass). In this particular case, however, 
almost all the galaxies for these high mass haloes will be 
brightest cluster galaxies, for which an al most flat colour- 
lumin osity relation is expected (see e.g. lEisenstein et al.l 
I2OOII) . in agreement with the model results. The scaling at 
this high-mass end goes as L oc M'^'^^ . This is similar to 
that obtained in other stud ies using the condition al luminos- 
ity function formalism, like lCoorav fc Milosavlie vic (2005a), 
who find that the central galaxy luminosity scales as ~ A/"'^ 
for high halo mass, while the total luminosity of the galax- 
ies in t he halo scale s as_~ M'^'^^ in the K-band; using 2dF 
results, lYang et al.l (|200^ also find that the luminosity of 
the central galaxy scales as Lc ~ M°'^^ . A fit to the obse r- 
vational K-band data gives L oc M°-^^ llLin fc Mohill2004ll . 

For low mass haloes and subhaloes, there are large dif- 
ferences between the results obtained by using the differ- 
ent luminosity functions. These are mostly a product of the 
characteristic luminosity, which determines where the break 
in the mass luminosity relation occurs, and the faint end 
slope, which determines the slope in the relation for low 
mass. These results illustrate the importance of an accurate 
determination of the faint end slope to obtaining a good 
mass luminosity relation. Nonetheless, looking at the results 
for the group luminosity it is possible to see that these dif- 
ferences have only a small effect when calculating the total 
luminosity associated with a high mass halo. The different 
results are quite close, particularly in terms of scaling, if 




Figure 6. Comparison of the galaxy (upper panel) and 
group (lower panel) luminosities obtained when using the non- 
parametric model described in section 4 with the different lumi- 
nosity functions of table 1. The luminosities have been scaled by 
the characteristic luminosity of the luminosity function used. 



less so than in the individual galaxy case, and the differ- 
ences seem mostly due to the characteristic luminosity i,. 
There is in fact some observational evidence that the scaling 
of cluster luminosi ty with mass is indep endent of photomet- 
ric band (see e.g.. iPopesso et al]l2004^ . Since the brightest 
galaxy luminosity gives a small relative contribution for high 
mass systems (cf flgurelKJ, it is possible to conclude that the 
main contribution to the luminosity should come from rel- 
atively massive subhaloes with luminosities roughly above 
O.IL*, below which the differences due to the low end slope 
become quite significant, mostly due to different faint end 
slopes in the luminosity functions. More importantly, this 
means that the mass-total luminosity relation obtained by 
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Band 0, (lO-^h^Mpc-^) a M, - Slogfc Ref. 



K 1.16 ±0.10 -1.09 ±0.06 -23.39 ±0.05 2MASS (Kochanek et al. 20011 

bj 1.61 ±0.08 -1.21 ±0.03 -19.66 ±0.07 2dFGRS fNorberg e t al. 2002) 

o-lg 2.18 ±0.08 -0.89 ±0.03 -19.39 ±0.02 SPSS (Blanton et al. 2003'! 

°-^r 1.49 ±0.04 -1.05 ±0.01 -20.44 ± 0.01 SPSS (Blanton et al. 20031 



Table 1. Schechter function parameters of the luminosity functions used. 



the non-parametric model depends mostly on the more well 
determined of the luminosity function parameters, L,. 



5.2 Cosmological parameters 

Although the agreement between our results and the large 
observational data set of brightest cluster galaxy luminosi- 
ties and halo masses shown in figure |S] was fairly good, it is 
possible to see that our model tends to slightly underpredict 
the expected luminosities. At such high masses, the subhalo 
contribution to the total host number and therefore to the 
mass-luminosity relation is neglegible (see figure|lj, and, as 
we have shown in the previous section, the results are also 
largely independent of the luminosity function. Therefore, 
the relation between halo mass and hosted brightest galaxy 
luminosity depends essentially on the halo mass function, 
and through it on the cosmology considered. This makes it 
important to look at the background cosmology used and 
how it affects the results. 

An example of this is shown in figure |7| where we also 
show the observational data for comparison. In general, it 
is possible to conclude that, comparing with the observa- 
tions, our model tends to prefer lower values of fim and as- 
The technical reason for this is quite straightfoward: lower 
values result in fewer high mass haloes, which means that 
when comparing to observed galaxy numbers it results in 
higher luminosity galaxies being associated with the same 
halo mass. For this reason, we have adopted in our calcula- 
tions a cosmological model with i}m ~ 0.25 and erg — 0.8, 
rather than the more standard Q,rn ~ 0.3 and erg — 0.9. 
This combination of values lies roughly on the lower edge of 
th e concordance region from WMAP+SDSS measurements 
of iTeemark et al ] i2004r ), and n ear the center of that from 
CMB+2dF dSanchez et al]|2005^ (see figure OJ. 

In fact, it is possible to obtain an even better agree- 
ment with the plotted observational values if we lower erg 
further. There a re two main rea s ons fo r not adopting such a 
model: first, the lTeemark et all i2004l results represent the 
best current measured estimate for the combination of the 
two; second, and more importantly, there are further addi- 
tional sources of error which are probably more significant 
than what is now a very accurate determination of the cos- 
mological parameters. The first of these is the possibility 
that the main assumptions of the model are wrong. Mainly, 
that is that we cannot treat the mass luminosity relation 
as one-to-one and monotonic. While this is certainly true 
in specific cases and we expect there to be significant scat- 
ter about the relation, it seems a fair assumption to make 
for the average relation as we are determining here; going 
beyond this assumption is outside the scope of this model. 
Second, and more importantly in the context of the present 
paper, there is the possibility that the model is giving in- 
correct results simply because the number functions we are 




a^=0.25, ag=0.8 

£1^=0.3, ag=0.9 

Lin & Mohr 2004 » 




Figure 7. Perived mass-luminosity relation for galaxies (upper 
panel) and groups (lower panel) as a function of the cosmological 
model used, for Qm = 0.25, erg = 0.8 and Qm, = 0.3, erg = 0.9. 

using as a basis for the comparison are mismatched, that is 
w(m)obsGrvcd 7^ n{m)modci, where the observational term ap- 
plies to the data points we are comparing our results with. 
Since the method is based on counting numbers for each 
of these, a difference between them would cause it to give 
incorrect results when comparing to observations. As we dis- 
cussed above, the cosmological model used can change the 



The Non-Parametric Model 9 




Figure 8. Group luminosity obtained by varying the low mass 
slope of the subhalo mass function, equation O between a = 1.5 
and a = 1.9 (our base model). 

model part of the equation. Far more likely, though, is that 
we are using an incorrect estimate for the observed mass that 
is skewing the observed term. We discuss this issue further 
in the next section. 

It is also worthwhile noting that, once the cosmology is 
fixed, the main determinant for the form of the halo mass- 
total luminosity relation is the slope of the subhalo mass 
function. This assumes that its normalization and cutoff are 
predefined, as discussed in section 3.2. Figure |H| shows the 
result of varying this slope. Generally, a flatter slope for the 
subhalo mass function results in a flatter slope at high halo 
mass for the total luminosity as well. Since a fit to the ob- 
servational data gives a slope of 0.72, while we obtain 0.88, 
using a flatter slope for the subhalo mass functi on than 1.9 
(quite possible in light of simulation results; see lShaw et alJ 
[2005) would likely give a better agreement with observa- 
tions, depending on possible corrections to their mass esti- 
mates, which may not be uniform. It should be noted, how- 
ever, that this in fact applies to the original subhalo mass 
function slope. How it compares with the present slope de- 
pends on the mass loss factor. With the one we are using in 
this paper, the slopes are mostly equal (see section 3). 

5.3 Observed mass determination 

As noted, one of the possibilities that may explain the differ- 
ence between our model results and the set of observational 
cluster data as presented in figure |S] is a potential misesti- 
mation of the corresponding halo mass. First of all, there 
is the problem of how to ext rapolate to the virial mass. In 
iLin. Mohr fc Stanford! i2004) . the authors determine M500, 
the mass enclosed in a radius where the average density is 
500 the critical, from X-ray temperature measurements. Be- 
sides the actual observational error in the temperature, and 
possible errors in the M500 — Tx relation used (which are as- 
sumed to be accounted for in the observational error bars). 




Figure 9. Halo mass function for our base model, together with 
the corresponding mass function for clusters as derived from the 
observed luminosity function and cluster halo mass-central galaxy 
luminosity relation (se e text for de t ails); the central line comes 
from the fit found by iLin fc Moliil i2004f) . the lines above and 
below represent the Icr range from the fit. Also shown is the curve 
for the flrn = 0.3, trg = 0.9 halo mass function, for comparison. 



there is the problem of how to extra polate M500 to the halo 
virial mass. Following what is done bv lLin. Mohr fc StanfordI 
i2004l) . we assume an NFW halo with a concentration of 5 
to calculate the virial mass, except that we take the virial 
radius to be at rioo rather than r2oo. A different density 
profile or simply concentration would result in a different 
measure for mass. Changing the concentration to a value in 
the range of between 4 and 7, for example, could result in as 
much as a 10% shift upwards or downwards in the estimated 
mass. 

The non-parametric model cannot work if the halo mass 
fucntion for the assumed cosmology does not match observa- 
tions, but there is a relatively straightforward way of check- 
ing how well these mass estimates agree with the results 
from our model. As mentioned previously, that is to check 
whether the predicted galaxy number at a given halo mass 
from the observed values match the number in the model, 
w(m)obsorvcd = n{Tn)^odoi, where nmodoi(m) in the high mass 
range we are interested in here is given by the halo mass 
function. As noted, these two must match if the comparison 
between the model and the data is to be meaningful, since 
the model is based on counting numbers. In order to deter- 
mine n(m.)observed, we combiue the luminosity function with 
the fit to the observed halo mass-br i ghtest cluster galaxy 
luminosity relation from lLin fc MohH (|2004) . The resulting 
mass function, together with the possible range obtained 
from the errors in the fit parameters, is shown in figure |n| 
From this we can see that, while the halo mass function 
we are using does not quite match the mass function corre- 
sponding to the observations, it is well within the range al- 
lowed by the fit errors. If we assume that the results from the 
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Figure 10. Mass to fej-band light ratio of galaxy systems (from 
isolated galaxies, to groups and clusters with growing mass). Data 
points are taken from the corrected mass to lig ht ratio calculate d 
for the 2PIGG catalogue of the 2dF survey bv lEke et all <200d) . 
These include both the directly measured values for luminosities 
above L > 10^'^h~^ Lq, and the values inferred by comparing 
the derived group luminosity function with a theoretical mass 
function, for luminosities below this threshold . The latter wer e 
extracted from the shaded region in figure 15 of lEke et al]i200d) . 
where the point represents the value of their model curve at that 
luminosity. The errorbars include uncertainties related to different 
possible values of erg . For comparison, we also show a li ne for the 
result of model D of lvan den Bosch. Yang &: Mol l2003l) . 



model are essentially correct, we can then use these curves to 
calculate a further correction to the estimated mass. In fact, 
this corresponds to finding the transformation that matches 
the fit to the observed data to our result for the average 
mass-luminosity relation. 



6 ADDITIONAL CHECKS 

In this section we calculate an additional four results which 
can be derived directly from our non-parametric model, and 
which can be used as an additional check on the success of 
the model; the mass to light ratio, the occupation number, 
the group luminosity function and the multiplicity function. 



6.1 Mass to light ratio 

The discussed trends of the mass- luminosity relation are fur- 
ther visible in figure 1101 which shows the &j-band mass to 
light ratio of the entire system (that is, a group or cluster if 
massive enough; this is basically a different way of plotting 
the result of figure EJ. 

As mentioned above, at the high mass end the clus- 
ter luminosity is almost directly proportional to halo mass 
(roughly, L oc A/"'****). This means that the resulting mass 
to light ratio will be almost constant, as can be seen in 



the figure, rising only very slowly with halo mass, which 
matches well with pr evious results for the mass to light 
ratio of clusters (e.g., 'Bahc all et aIll200Ct iKochanek et alJ 
,2003 ; Eko ct al. 2004. .2005.,) . In this case, though, the val- 
ues we obtain for the mass to light ratio seem to be 
slightly smaller than the observational results at the bright 
end, otherwise they seem in good agreement. This is re- 
flected in the values for the cluster mass to light rat io. 
The value derived bv iFuku gita^ Hogan fc Peebles! l|l99dl is 
450 ± lOO/i(M/L)0, while lEke et al.l ll2004ll obtain an aver- 
age value of 466 ± 2Qh{M / L)q; we obtain a slightly lower 
cluster mass to light ratio of approximately 425h{M / L)q 
for a lO^^h-^Me halo. Since we obta ined a good agreement 
with t he cluster luminosity results of lLin. Mohr fc StanfordI 
i2004h and, as discussed above, the mass luminosity for sin- 
gle galaxies is well constrained in our model, the most likely 
origin for the discrepancy we find here is likely to be the sub- 
halo distribution. Lowering the value of the subhalo mass 
function slope slightly results in a better agreement with 
the observational results at both the low and high luminos- 
ity ends; however, the agreement is poorer at intermediate 
luminosities, while the slope of the derived mass to light 
ratio becomes steeper (cf. results in figure |HJ. 

We can also comp are our results with the conditiona l 
luminosity function of lYang. Mo fc van den BoschI J2003D . 
Plotted in figure 1101 is the result of the average halo 
mass to light ratio fit with the p arameters of model D of 
Ivan den Bosch. Yang fc Mol J2003t) . There is a good agree- 
ment at the low end, while again we find some disagreement 
at the high end. Note that this disagreement is most likely 
due to the fact that we are adopting a different cosmol- 
ogy than these authors, with slightly lower values for Qm 
and (Tg. Nevertheless, the overall shape is quite similar; the 
small difference is most likely due to the factor that these 
authors are fitting the mass to light ratio to a double power 
law, with a sharp transition to a constant at high mass, 
while we obtain a smooth function. The minimum also oc- 
curs at a slightly higher luminosity. We are actually in better 
agreement with the observed values at an intermediate range 
(around 10^" — IO^^/i^'^Lq), while the opposite happens at 
the bright end. between both results is considerably better; 
as noted, using a flatter slope for the subhalo mass function 
may help with the problem. 

Figure 1101 also illustrates quite well the various mass 
(or, as presented in the figure , luminosity ) scales which de- 
termine galaxy properties (cf. (Dete^HOQa)- There is a mini- 
mum in the mass to light ratio at a mass scale of around 
3 — 4 X 10^^/i~^Mq. This is a good match for the scale 
which represents a shift in the characteristic galaxy popula- 
tion corresponding to a change in the s ources of gas accre- 
tion a nd star formation suppression (e.g.. lDekel fc BirnboimI 
l2004h . Also, as we have already shown, the steep slope 
in the mass luminosity relation means that haloes below 
roughly a few times lO®/i~^M0 will essentially be dark; 
again, this matches well with the mass belo w which photo- 
ionization is expected to suppress gas infall jlBabul fc Ree3 
119921: iThoul fc Weinberelll996l: iDekel fc Wocll2003H . At the 
high end, there is a noticeable break in the mass to light ra- 
tio at a corresponding mass of around 10^^/i"^Mq. As can 
be seen in figure |51 this marks the scale where we go over 
from isolated galaxies with at most very small satellites to 
actual groups, and once more is in good agreement with 
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what is observed and expected from a theoretical point of 
view, which predicts an upper bound for cooling in a dy- 
namical time at this mass scale. This is also the mass scale 
at which merginR is expected to become inefficient (e.g., 
ICoorav fc Milosavlievi3l2005all . 



6.2 Occupation number 

With the subhalo mass function given by equation it is 
straightfoward to calculate the occupation number (that is, 
the number of subhaloes in a parent halo of given mass), as 
a function of halo mass M. This is simply given by: 



1000 



Ns{M) 



N{m\M)dm. 



(10) 



The only complication is that it is necessary to specify a 
minimum mass for the subhaloes, nimin', otherwise, the in- 
tegral is divergent. Since we equate subhaloes with satel- 
lite galaxies, we can associate this minimum mass with a 
minimum luminosity of these galaxies, and then Ns{M) + 1 
will give us the total number of galaxies in a halo of mass 
M, with luminosity greater than the minimum we are con- 
sidering; this is one of the key ingredients of the HOD 
models, (e.g. iBerlind fc Wein berg 2002; Zhen g et al J 120041: 
lYaneet al.l2005^ . These models usually take the function in 
mot to be the average value of the number of satellite galax- 
ies, with the actual number a Poisson distribution around 
this average; this seems to be in good agreement with simula- 
tion results (Kravtsov ct al. 2004; Zheng ct al. 2004). Also, 
the occupation number of the central galaxy is usually taken 
to be a step function, being or 1 depending on whether 
the halo mass is greater than a certain minimum; the case 
in our model is also similar, as we consider that a halo hosts 
a galaxy with luminosity above a certain threshold if its 
mass is greater than that corresponding to this luminosity 
from the mass-luminosity relation. Our results for the oc- 
cupation number are shown in figure [TTI where we take the 
minimum mass to be the necessary to host a galaxy with 
Mk = —21. The points plotte d are t aken from the cluster 
data of iLin. Mohr fc StanfordI i2004l) . where as before we 
have included a correction to the mass. 

Our result is in quite good agreement with the plot- 
ted data. At high halo mass (approximately above ~ 
lQ^^h~^ Mq), our result for the occupation number scales 
roughly as A'' oc A/"'^. Qualitatively, it also compares well 
with those obtained from simulations, analytical models and 
observations by a variety of authors. The first thing to note 
is that, unlike the result in paper I, the occupation number 
is now no longer a function of M/rrimin alone: there is an 
extra dependence on the halo mass M. Physically, this is un- 
derstandable as haloes having M/M^ > 1 are in the process 
of growing and merging whereas those having M/M^ < 1 
are decreasing in number density as they merge into larger 
systems. This is a consequence of the fact that the the sub- 
halo mass function, given in equation Q, cannot be written 
as a function of solely m/M: its normalization has terms 
which depend only on M. While a direct comparison of 
the numbers is complicated by the fact that they depend 
on what mass threshold is being considered, the slope of 
the occupation number for high mass haloes is roughly the 
same as that found in other studies. This matches well with 
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Figure 11. Occupation number: average number of galaxies pre- 
dicted for a parent halo of mass M, given by the sum of the 
number of subhaloes with galaxies above a given luminosity plus 
one for the central galaxy. The luminosity threshold adopted was 
M]f < —21, which corresponds to the minimum used for the 
observed data points plotted, taken from llTn. Mohr fc StanfordI 
120041 . 



the an alytically derived halo mass function of lOguri fc Led 
liooi), who find a slope of 0.9 at high M/m « lO'^, and 
closer to 1 at lower M/m « 10 ^, and also with the r esults 
from simulations: for example, [Kravtsov e^^ ^2004^ find 
a slope very close to 1, while lzh^n^^^£jTl200^ . in fitting 
HOD models to semi-analytic models of galaxy formation 
and smoothed particle hydrodynamic simulations, find val- 
ues for the slope between 0.97 and 1.24 (depending on the 
baryonic mass threshold). Results of obse rvational studies 
also seem to agree on a slope of about 1: iKochanek et alj 
(2003) obtain a relation for t he number of galax:ies with 
L > L, scaling as M}j^ , while IXbaz aiian et alJ ('2005') fit a 
HOD model, together with the cosmological parameters, to 
the projected correlation function of a volume limited sub- 
sample of the Sloan Digital Sky Survey (SDSS), together 
with CMB results, and find a good agreement with models 
where the slope is fixed at a value of 1; when they leave this 
as a free parameter, they obtain a result of 0.8, but without 
a significant improvement in the quality of the fit. 



6.3 Group luminosity function 

Also of interest when studying clusters is the group lumi- 
nosity function, (f)g{Lg). This is the analogous of the galaxy 
luminosity function, but is calculated for groups and clus- 
ters, and gives the number density of these objects in a given 
luminosity range. We can obtain this in our model by trans- 
forming the mass dependent halo mass function by using the 
relation between halo mass and total group luminosity: 



4)g{Lg)dLg 



nh{M{Lg))—dLg 



(11) 
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Figure 12. Group luminosity function. The solid red line rep- 
resents the result of the non-parametric model using equa- 
tion illl : the two other lines are observational group lumi- 
nosity functio ns extracted from gala xy catalogues, respectively 
the AGS I M oore. Frenk fc Whit3[l99a') (dotted) and the VSLF 
l lMarinoni. HudsorrirGiuricinll2002l) (dashed) . 

Here, we are implicitly assuming that a group/cluster is 
made up of two components: a central galaxy hosted in the 
parent halo itself, and its satellites, hosted by the subhaloes. 
Thus the total luminosity is the sum of these two contribu- 
tions: 

Lg{M)=L(M)+l L{m)N{m\M)dm. (12) 
Jo 

This group luminosity is the same as was presented in fig- 
ure|S]above. Our result for the group luminosity function is 
shown in figure The determination of the observational 
function is usually done from galaxy catalogues by building 
groups of gravitationally bound galaxies. In the figure, we 
show the luminosity function fits of two different such func- 
tions : the AGS for the CfA survey jMoore. Frenk fc Whit3 
Il993^ and the VSLF iMarinoni. Hudson fc GiuricinI l2002l) 
group luminosity functions. We also show data for the 
group luminosity function as measured from 2dF results of 
lEke et aP l|20o5. We have included corrections to the AGS 
results by using bj — -Bzwicky — 0.05 and making the VSLF 
results 0.55 magnitudes fainter to compensate for the differ- 
ence between the B and bj band s and internal absorption 
(see comments in lEke et al.ll20()5l) . We obtain a very good 
agreement with the observational results. This agreement is 
markedly better than the result in paper I; this is essen- 
tially a result of the higher group luminosity we obtain for 
each halo, which in turn is caused by an increase in the 
subhalo numbers at a given luminosity due to the more self- 
consistent way in which we are treating mass loss. 

6.4 Multiplicity function 

It is also possible to derive the multiplicity function, the 
number density of group/clusters as a function of their rich- 
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Figure 13. Multiplicity function derived from the non- 
parametric model (solid line), based on combining the halo mass 
function with the total halo occupation number, that is, the total 
number of galaxies, above a given luminosity, present in a parent 
halo. The minimum mass for the halo or any of its subhaloes to be 
counted corresponds to a magnitude of Mb = —19.4. The differ- 
ent po ints are taken from the results derived bv lPeacock fc SmithI 
120001) from the CfA survey, while the two additional lines are es- 
timates based on the group luminosity functions shown in figure 
1121 where the magnitude limit for all of these is also Mg = —19.4. 



ness, by a process similar to the one described above for the 
group luminosity function. To do this, we replace the group 
luminosity Lg in equation IIH with the number of galaxies 
present in the halo. This is given by the total occupation 
number calculated above and shown in figure 1111 that is, 
the number of luminous subhaloes plus one. The only other 
thing that is necessary to take into account is the limit to 
the galaxy luminosity we want to consider. This will then 
give a lower limit to the mass of the haloes and their sub- 
haloes to put into the occupation number calculation. Our 
results for the multiplicity function are shown in figure ITHl 
We have taken as the lower luminosity limit Mb = —19.4, 
in order to match the observational results also shown. Due 
to the way in which it is built, there is a sharp upturn at 
TV = 1 due to the haloes which are massive enough to con- 
tain a galaxy but not to have subhaloes big enough to host 
a galaxy themselves. This is in part a consequence of mod- 
elling the occupation number of the parent halo as a step 
function, as was done above; a more realistic model would 
have a smoother transtition (see e.g. Zheng et al. 2004 for an 
example of how this is done in the context of HOD models) , 
which would help to atenuate this. 

Figure 1131 also shows some observ ational data. The 
points were taken from the analysis of iPeacock fc SmithI 
^2000), while the two lines were constructed from the two 
group luminosity functions shown in figure [m To do this, 
we assume that the luminosity function of the galaxies in 
each group has the same shape as the general galaxy lu- 
minosity function of the survey, a Schechter function with 
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characteristic luminosity L« and faint end slope a, but with 
the normalization n determined by the group luminosity. By 
integrating this function, the group luminosity Lg can then 
be calculated as Lg = nL*r(2 — q); this is then related to the 
number of gala^des in the group above a certain threshold 
luminosity Lmin by 

^ = Tr^^^^0^^ (13) 

where we used the expression above to replace n. Finally we 
can use this relation between group luminosity and mem- 
ber numbers in the group luminosity function to obtain the 
multiplicity function, using n{N) = <j>g{Lg)dLg/dN , where 
n{N) and (f)g[Lg) are the multiplicity and group luminos- 
ity functions, respectively. The galaxy numbers for both 
the points and the results derived from the group luminos- 
ity functions were then multiplied by a factor of 0.66 to 
take into account the difference in radius between friends 
of fr iends estimates and th e usual definition of virial radius 
(see iKochanek et al.ll2003^ . 

We obtain a good agreement with the observational 
points, better than was the case in paper I. This is not sur- 
prising given the results we presented above for the group 
luminosity function, and the root cause is the same: the 
more self-consistent treatment of mass loss has resulted in 
an increase in subhalo numbers at a given mass. On the 
other hand, at high member count, our result would seem 
to be above the multiplicity functions calculated from the 
group luminosity function. This is most likely due to the 
fact that we assumed above that we can use for the lumi- 
nosity function of the group galaxies the same shape as the 
background luminosity function; while this should be true 
for moderately dense regions (corresponding to low A'^), it 
is not a very good approximation for v ery dense environ- 
ments like clusters, where A'' is high fsee lCroton et al']l2005l 
for observational results). In fact, a higher L», expected in 
clusters with high A'', gives a higher value for the multiplic- 
ity function, therefore bringing the two curves into better 
agreement with our results. 



7 CONCLUSION 

In this paper we have presented a non-parametric model to 
relate the luminosity of a galaxy to the mass of the halo or 
subhalo hosting it. We start by assuming that this relation 
is one to one and monotonic, and then compare their num- 
ber as given by their statistical distributions, obtained from 
observations and simulations, respectively. This gives us the 
average relation between the luminosity of a galaxy and the 
mass of the dark matter halo or subhalo which hosts it. We 
can also determine the total group luminosity as a function 
of halo mass by integrating over the luminosity of all the 
subhaloes in a given parent. 

We have argued that, to maintain the assumption of 
monotonicity, when accounting for the subhaloes it is in fact 
necessary to consider not their present mass, but rather their 
original mass just prior to accretion into their parent halo. 
In order to allow for this, we have introduced a simplified 
prescription for the average mass loss factor that we can 
then apply to subhalo mass functions measured from sim- 
ulations to regain the original subhalo mass function. We 



have noted that, however, there are two strong constraints 
on this function: first, no subhalo could originally have more 
than half the present mass of the parent halo, otherwise it 
would be, by definition, the parent halo itself; second, it 
is to be expected that, since the parent is built up by ac- 
creting and stripping mass off the subhaloes, that the total 
mass in the original subhaloes equals the present mass of 
the parent halo. If the mass loss factor is fairly regular (or 
more precisely, its logarithmic derivative is always signifi- 
cantly lower than 1, which is mostly the case with the one 
we presented), and assuming a present day Schechter type 
subhalo mass function (w hich seems to be obtained from 
simulation results; see e.g., IShaw et all2005ll . these two con- 
straints mean that the only free parameter in the original 
subhalo mass function would be the low mass slope, which in 
these circumstances would be close to the present day one. 

We find that, for high mass haloes, the mass- luminosity 
relation appears to be mostly independent of the luminosity 
function being used, when the luminosity is scaled to the 
characteristic luminosity L,. The same is true, to a rather 
lesser extent, of the mass-group luminosity relation, which 
has an increased dependence on L, . At the low mass end, the 
break in the relation is associated with L», while the slope 
depends significantly on the faint end slope of the luminosity 
function. 

Overall, our results are a good match to observations 
and results of other theoretical models. We find, however, 
that for high mass haloes, our results seem to slightly un- 
derestimate the luminosity of both central galaxies and 
clusters when compared to the observational resu lts of 
iLin. Mohr fc Stanford! (jgOOi) and lLin fc Mohrl 12004). This 
small discrepancy between our results and observational 
data simultaneously raises some concern. The reason is that, 
as we have discussed, the mass luminosity relation in high 
mass haloes is practically independent of both the actual lu- 
minosity function used and of the subhalo population. This 
means that our result depends only on the halo mass func- 
tion, and through it, on the cosmological model. We have 
shown that, when comparing our results with observations, 
our model seems to prefer lower values of Q.m a nd erg ; within 
the concordance region of lTeemark et alJ i2004ft . best results 
are obtained near the lower boundary, and consequently we 
have used Q,m ~ 0.25 and as = 0.8 to construct our base 
model. While to change our result would necessitate chang- 
ing the cosmology, or more seriously, the central assumptions 
underlying the model, it is far more likely that the discrep- 
ancy is actually a product of a misestimation of the observed 
halo mass. In fact, there is some uncertainty regarding the 
values of the mass cited in these studies, specifically in ex- 
trapolating them to the virial radius. A reasonable change in 
the concentration used to make this extrapolation may ac- 
count for as much as a 10% shift in the estimated observed 
mass. Further, using the fit to the observed mass- luminosity 
data and the luminosity function, we have compared the ex- 
pected number of galaxies in haloes of a given mass to the 
one we are using. We find that, while the agreement is not 
perfect, the mass function we used is well within the error 
range of the observed one. Further, assuming that our model 
is correct, it is possible to view the required transformation 
to the estimated mass to match the two curves as an addi- 
tional correction to it. 

The situation with the total luminosity on cluster scales 
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is relatively similar, though here the problem might be more 
one of shape rather than normalization. The slight discrep- 
ancy between our results and the observational data may be 
improved if the halo mass has been underestimated obser- 
vationally or by tweaking with the cosmological parameters. 
There is however an additional factor, in that the subhalo 
population contributes significantly to the total luminosity. 
This is represented by the original subhalo mass function, 
and we have argued that its normalization and cutoff mass 
should be considered fixed, the former to give the total 
present parent halo mass, the latter to avoid having sub- 
haloes which originally had more than half the parent mass. 
If we assume that this function has a Schechter shape (as 
would be the case with a present subhalo mass function with 
a Schechter shape and a regular mass loss factor), then the 
only free parameter left is the low mass slope. We have also 
argued, by comparing the results obtained using different 
luminosity functions, that the more important luminosity 
function parameter that determines the result at this high 
mass end is the characteristic luminosity L, , while the faint 
end slope (which could be considered the most uncertain of 
the parameters) does not seem to be as big a factor. There- 
fore, if we take the subhalo mass-galaxy luminosity relation 
(and consequently the cosmology) to be fixed, the remain- 
ing free parameter in the model is the subhalo mass function 
low mass end slope. 

On the one hand, this means that we cannot much vary 
the total number of satellite galaxies in this framework. On 
the other hand, our result for the group luminosity seems to 
be a better match to the observed values than that of the 
galaxy luminosity; if we were to increase the latter, we might 
end up obtaining too high a luminosity in groups, assuming 
the slope of the subhalo mass function is kept fixed. This po- 
tential problem may however be viewed in light of the recent 
results of Coorav & Ccn ( 2005) , who claim that when using 
a similar prescription to model the satellite galaxy luminos- 
ity function using the subhalo mass function, it is necessary 
to introduce an efficiency function lowering the number of 
luminous subhaloes otherwise predicted to match the two. 
There is also the possibility, which we have not considered, 
that the mass-luminosity relation is different for haloes and 
subhaloes; using a different relation for the latter might ex- 
plain the difference to the observational results. Since gas 
accretion and mergers stop in subhaloes once they are ac- 
creted, we are led to the conclusion that, if this were to be 
the case, the subhaloes would actually be less luminous than 
what we are considering. Additionally, there is the effect of 
the slope of the subhalo mass function. Using a flatter slope 
than the 1.9 value we have used for our base model may 
help with this situation, since it would flatten the slope of 
the group luminosity calculated without changing the high 
mass galaxy luminosity. This may be bcnificial, since the fit 
to the observed clusters fromlLhi. Mohr & Stanford (20o3) 
goes as L oc M"'^'^, while we get a slightly steeper 0.88. If 
we do accept the whole framework of the model, though, 
this raises an intriguing possibility: the only free parameter 
we have is the low mass end slope of the subhalo mass func- 
tion, and so we may be able to get a completely independent 
confirmation of its value by fitting the model results to the 
observational data. 

Quantitavely, we find that central galaxy luminosity 
scales with halo mass as Li ~ M°-2* for high mass haloes. 



fairly independently of waveband when the luminosity is 
scaled by the appropriate characteristic luminosity L,, and 
also of the form of the subhalo mass function. The total 
group luminosity scales with halo mass as Ltot ~ A/" ***, 
also fairly independently of waveband when appropriately 
scaled; a flatter subhalo mass function low mass slope results 
in a flatter dependence at the high mass end. This implies 
that the halo mass-to-light ratio is almost flat at high mass 
(and luminosity), and we obtain a value of 425h{Mi, j / L)q 
for a lO^^h'^AlQ halo. For low mass haloes, the resulting 
mass luminosity relation is dependent on what waveband 
we are considering, but scaled luminosity goes as L ~ 
with a roughly between 4 and 4.5 for haloes or subhaloes in 
the mass range 10^ Mq to 10^° Mq; for example, in 
the K-band used for our base model, a ~ 4.5. We also flnd 
that the occupation number, which can be derived almost 
directly from the subhalo mass function, scales as M'^'^ . 

Finally, we should also make a note on the applicabil- 
ity of the relation we obtained. First, we should stress that 
this is an average relation. As can be seen from the plotted 
observational data, we expect a rather large scatter around 
it. While we feel that obtaining this average relation is quite 
an important flrst step and by itself already allows a range 
of applications, it is important to obtain a model for the 
scatter if it is to be applied, for example, to build mock cat- 
alogues from simulation results. In this direction, there has 
already been some work on applying the base framework 
that we have developed further here to a context of a con- 
ditional luminosity function, including potential prescrip- 
tions for scatter (see work by Coora y and collaborators, e.g. 
ICooravllioOSbl: ICoorav fc Cenll2005^ . Nonetheless, we feel it 
is relevant to have a good analysis of the basic framework, 
especially since the overall simplicity of the model makes it 
conceptually clear and pedagogical, while at the same time 
allowing a good comprehension of the factors influencing it. 
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